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ABSTRACT 

We present a survey of 2 ; ^ 3 Lyct emitters (LAEs) within the fields of the VLT LBG 
Redshift Survey. The data encompasses 5 independent survey fields co-spatial with 
spectroscopic LBG data and covering a larger total area than previously analysed for 
LAE number counts and clustering. This affords an improved analysis over previous 
work by minimising the effects of cosmic variance and allowing the cross-clustering 
analysis of LAEs and LBGs. Our photometric sample consists of ~ 600 LAE candi¬ 
dates, over an area of 1.07 deg^, with equivalent widths of > 65 A and a flux limit of 
^2x10“^^ erg cm“^ s“^. Erom spectroscopic follow-up, we measured a success rate 
of 78 ± 18%. We find the R-band continuum luminosity function to be ^ 10x lower 
than the luminosity function of LBGs at this redshift, consistent with previous studies. 
Exploiting the large area of the survey, we estimate the LAE auto-correlation function 
and find a clustering length of ro = 2.86 =b 0.33 h~^ Mpc, low compared to the 2 ; ^ 3 
LBG population, but somewhat higher than previous LAE measurements. This corre¬ 
sponds to a median halo mass of Mdm = 10^ h~^ M©. We present an analysis 
of clustering length versus continuum magnitude and find that the measurements for 
LAEs and LBGs are consistent at faint magnitudes. Our combined dataset of LAEs 
and LBGs allows us to measure, for the first time, the LBG-LAE cross-correlation, 
finding a clustering length of ro = 3.29 zb 0.57 h~^ Mpc and a LAE halo mass of 
10 ii i^O-4 ^-1 Mq. Overall, we conclude that LAEs inhabit primarily low mass halos, 
but form a relatively small proportion of the galaxy population found in such halos. 

Key words: (cosmology:) large-scale structure of Universe, galaxies: evolution, galax¬ 
ies: high-redshift, galaxies: luminosity function, mass function 


1 INTRODUCTION 

Two key methods for efficiently selecting high redshift galax¬ 
ies are via their Lyman break feature (Lyman Break Galax¬ 
ies - LBGs) and via their Lyman a emission (Lyman-o 
emitters - LAEs). Selecting high redshift galaxies through 
their strong emission in the Lyo feature using the narrow- 
band imaging method has come to be a very effec¬ 
tive technique to isolate high redshift galaxies. There are 
many observations that have been made to uncover galax¬ 
ies with strong Lya emission at various redshifts (e.g. 
Hu & McMahon 1996; Hu et al. 1998; Gowie & Hu 1998; 
Steidel et al. 1996; Ouchi et al. 2003; Hayashino et al. 2004; 
Gawiser et al. 2007; Gronwall et al. 2007; Ouchi et al. 2008; 
Rauch et al. 2008; Ouchi et al. 2010; Bahados et al. 2013). 
At redshift 2 ; ~ 3, there are a number of photometric 
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and spectroscopic LAE samples (e.g. Steidel et al. 2000; 
Fynbo et al. 2003; Hayashino et al. 2004; Matsuda et al. 
2005; Venemans et al. 2007; Nilsson et al. 2007; 
Gawiser et al. 2007; Gronwall et al. 2007; Rauch et al. 2008; 
Ouchi et al. 2008; McLinden et al. 2014). Gronwall et al. 
(2007) surveyed the Extended Chandra Deep Field-South 
(ECDFS), identifying 160 photometrically selected 
2 ; = 3.1 LAEs over an area of 0.28 deg^ and to a Lya 
luminosity limit of LLya ^ 10^^ erg s“^. Their colour 
constraints correspond to a rest-frame equivalent width 
limit of EW > 20 A. By measuring the Lya and rest-frame 
UV continuum luminosity functions in combination with 
the equivalent width distribution, they surmised that LAEs 
contain a non-negligible amount of dust and are therefore 
not of a primordial origin (but that they do represent a 
young galaxy population). Ouchi et al. (2008) presented 
a survey of ~ 350 photometrically identified (with 41 
spectroscopically confirmed) LAEs over 1 deg^ in the 
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Subam/XMM-Newton Deep Survey (SXDS) field to a 
limiting Lya luminosity of L^ycx ^ 10^^ erg (and a 
rest-frame equivalent width limit of > 65 A). Combining 
their z = 3.1 data with higher redshift samples, they find 
that LAEs (at a given equivalent width and luminosity 
limit) are more common at earlier epochs and are a tracer of 
either low-extinct ion (when compared to LBGs) or younger 
stellar populations. 

Although LAEs in general have been widely studied, 
their clustering properties remain comparatively sparsely 
studied, due to the need for large survey fields. At 
z ^ 3, Hayashino et al. (2004), Gawiser et al. (2007) and 
Ouchi et al. (2010) have each presented measurements of the 
LAE angular correlation function. Hayashino et al. (2004) 
do not present a quantification of their clustering result, 
whilst Gawiser et al. (2007) and Ouchi et al. (2010) mea¬ 
sure clustering lengths of ro = 2.32^^^'^ h~^ Mpc and 

ro = 1.70lQ 4 g Mpc respectively. These two measure¬ 
ments are based on differing equivalent width constraints 
of EW > 20 A and EW > 65 A respectively, with the 
lower EW cut giving a larger clustering length (although the 
two measurements only differ by ~ 1 cr). At higher redshifts, 
LAE clustering measurements have been obtained at z ~ 4 
(Kovac et al. 2007; Ouchi et al. 2010), z ~ 5 (Ouchi et al. 
2003; Shimasaku et al. 2004; Shioya et al. 2009), and z ~ 
6 — 7 (Ouchi et al. 2010). In each case, these clustering mea¬ 
surements are based on single observational fields of at most 
1 deg^, whilst in some cases (e.g. Hayashino et al. 2004, at 
z = 3.1) they are actively targeted on high density regions. 
Indeed, cosmic variance can be a significant source of error 
in LAE clustering studies (Shimasaku et al. 2004; Erancke 
2009) and measurements based on single fields can prove 
unrepresentative of the population as a whole. The largest 
of these works thus far remains that of Ouchi et al. (2010), 
which presented LAE clustering results atz = 3.1, z = 3.7, 
z = 5.7, and z = 6.6. An apparent evolution in the LAE 
clustering was observed, with ro (and therefore the clus¬ 
tering bias) increasing with increasing redshifts, implying a 
constant host halo mass with redshift of ~ 10^^ /i“^Mpc. 

Here we present observations of LAEs in the deep imag¬ 
ing fields of the VLT LBG Redshift Survey (VLRS) aimed 
at studying the clustering properties of LAE galaxies. Our 
LAE sample extends over 5 independent observational fields 
each measuring ~ 0.5° x 0.5° in area, thus combatting the ef¬ 
fects of cosmic variance. These fields were selected based on 
the presence of bright quasars at z > 3 and so are not specif¬ 
ically targeted at galaxy overdensities as with some previous 
works (e.g. Hayashino et al. 2004). We target the z ^ 3 LAE 
population in order to overlap with our z ^ 3 spectroscopic 
LBG sample in these fields, allowing cross analysis between 
the two populations. The primary aim of the VLRS has been 
to probe the intergalactic medium (IGM) of z ^ 3 galaxies 
by observing galaxies close QSO sightlines (Grighton et al. 
2011; Tummuangpak et al. 2014). At present the survey con¬ 
sists of 9 fields in total, with spectroscopic observations of 
Lyman break galaxies (LBGs) taken with VLT VIMOS. The 
primary survey is described in detail by Bielby et al. (2011, 
2013). By adding narrow-band data to our survey data, we 
aim to enhance the scope of the survey, probing the volume 
for both LAEs and other sources of Lya photons (e.g. Lya 
blobs and Lya emission from and around QSOs). 

In this paper, we present the narrow band data and 
analysis of LAE galaxies within the VLRS. The data covers 
5 independent survey fields, giving the significant advan¬ 
tage of reducing the potential impact of cosmic variance on 
our results compared to past LAE surveys. We discuss our 
deep imaging data in these fields and the selection of LAEs 
at z ~ 3.1 using the NB497 narrow band filter in Sec. 2. 


We then report on spectroscopic follow up observations and 
our optimised selection criteria based on these in Sec. 3. 
In Sec. 4, we present narrow band number counts and the 
rest-frame UV continuum luminosity function for our LAE 
sample, whilst in Sec. 5 we present an analysis of the clus¬ 
tering of the LAEs. Einally we present a discussion of our 
results in Sec. 6 and our conclusions in Sec. 7. 

Throughout this work, we adopt a cosmology given by 
(Qrn, h, as) = (0.25, 0.75, 0.73, 0.83), where h is the Hub¬ 
ble constant in units of 100 km/s/Mpc. Distances and vol¬ 
umes are given in comoving coordinates unless otherwise 
stated. All magnitudes are given in the AB system unless 
otherwise stated. 


2 IMAGING DATA 

The VLRS consists of 9 deep fields (each measuring be¬ 
tween 0.25 deg^ and 1 deg^) centred on bright z > 3 QSOs 
and containing both deep broad band photometry and spec¬ 
troscopy (Bielby et al. 2011, 2013). Here we have added to 
these data by using the Subaru Suprime-Gam instrument 
(Miyazaki et al. 2002) to obtain deep narrow band (497 nm) 
imaging in 5 of the VLRS fields, which allows for identifica¬ 
tion of LAEs at z ~ 3.1 (details in Sec. 3.2). The observa¬ 
tions were taken 19th September 2009. 

Suprime-Gam is a mosaic GGD camera with ten GGDs, 
measuring 2048 x 4096 pixels, which covers a 34' x 27' field 
of view with a pixel scale of 0.20". We observed ~ 0.5° x 
0.5° fields centred on the bright quasars QSO B2359+068, 
LBQS 0301-0035, QSO J0124+0044, PKS 2126-158, and 
LBQS 2231—0015 (details are in Table 1). These fields were 
observed with the narrow band [OHI] NB497 filter (4977 A, 
EWHM 77 A; Hayashino et al. 2003). Individual expo¬ 
sure times were 1,200 s for each frame and the seeing was 
generally sub-arcsecond (~ 0.6" — 0.7"). Eor the purposes 
of photometric calibration, we observed the standard star 
LTT 9491 (R.A. 23:19:34.98, Dec. -17:05:29.8 J2000). 

In order to select LAE galaxies, we also require B 
and V or R deep broad band imaging, which are avail¬ 
able from our previous work in these fields. The full de¬ 
tails of the broad band data and its reduction are provided 
by Bouche & Lowenthal (2004); Bielby et al. (2011, 2013), 
whilst a summary of the broad band data is given in Table 2. 


2.1 Data Reduction 

To reduce the Suprime-Gam raw narrowband data, we 
used the pipeline software, sdfred (the Suprime-Gam 
Deep field REDuction package^, Ouchi et al. 2004a) which 
comprises iraf, SExtractor (Bertin V Arnouts 1996), 
and the mosaic-GGD data reduction software (Yagi et al. 
2002). The package includes bias subtraction, flat field¬ 
ing, distortion+atmospheric dispersion corrections, match¬ 
ing the PSE size, sky subtraction, masking vignetting caused 
by the auto-guider, masking bad pixels, image alignments 
and scaling, and mosaicing. The reduction procedure is 
briefly described below. 

The images were first bias subtracted using the me¬ 
dian value of the overscan region on each line of a given 
file (the Suprime-GAM GGDs are noted to have very little 
bias pattern, so only the subtraction of the overscan median 
is required) and the overscan regions were then trimmed 
from the images. The flat field was then created using a 
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Table 1. Details of narrowband imaging (NB497) from the Subaru/Suprime-Cam. 


Central QSO 

R.A. Dec. 

(J2000) 

QSO redshift 

^exp 

(S) 

Seeing 

(") 

Mzp 

la depth 

QSO B2359+068 

00:01:40.57 

+07:09:54.1 

3.234 

6000 

0.61 

31.94 

28.38 

QSO J0124+0044 

01:24:03.78 

+00:44:32.7 

3.807 

4800 

0.55 

31.93 

27.47 

LBQS 0301-0035 

03:03:41.05 

-00:23:21.9 

3.175 

7200 

0.60 

31.93 

28.57 

PKS 2126-158 

21:29:12.18 

-15:38:41.0 

3.268 

6000 

0.96 

31.89 

27.90 

QSO B2231-0015 

22:34:08.99 

+00:00:01.7 

3.027 

6000 

0.87 

31.93 

27.77 


Table 2. Details of broadband imaging observations from Bielby et al. (2011, 2013). 


Central QSO 

Band 

'Vexp 

(ks) 

Seeing 

(arcsec) 

la depth 

Instrument 

QSO B2359+068 

B 

7.2 

1.45 

28.23 

KPNO/MOSAIC 


R 

6.0 

1.15 

27.76 


QSO J0124+0044 

B 

2.8 

1.50 

27.99 

KPNO/MOSAIC 


V 

3.1 

1.40 

27.74 


LBQS 0301-0035 

B 

6.4 

1.28 

28.30 

KPNO/MOSAIC 


R 

4.8 

1.19 

27.66 


PKS 2126-158 

B 

7.8 

1.60 

28.44 

CTIO/MOSAIC2 


R 

6.4 

1.50 

28.27 


QSO B2231-0015 

B 

13.2 

1.01 

27.28 

WFCam (INT) 


R 

19.2 

1.01 

26.98 



total of 25 dithered object frames from across the 5 ob¬ 
served fields, with objects and regions vignetted by the auto- 
guider (AG) probe masked. With this flat field applied to 
the individual object frames, the astrometric distortion cor¬ 
rection was performed (correcting for the telescope optics 
and the differential atmospheric dispersion) using a 5th or¬ 
der polynomial transformation (Miyazaki et al. 2002). The 
point spread function (PSF) was then matched across the 
images by applying a gaussian smoothing kernel to individ¬ 
ual images. The sky background was then determined by 
interpolating over a mesh pattern. The image was divided 
into a grid of 64 x 64 pixel squares (corresponding to 12'^9 
X 12'^9). A bilinear interpolation was used to determine the 
global sky background from this grid and the result was 
subtracted from the individual science frames. 

An initial astrometry solution was calculated using SD- 
FRED, whereby the images were matched internally (i.e. us¬ 
ing one of the science frames as a reference) and matched 
to a reference stellar catalogue. The images from each field 
were then co-added into stacked final images. 

We applied the cosmic ray rejection using the rejected- 
mean algorithm, crreject, from IRAF. 

2.2 Object Detection and Photometry 

Extraction of sources from the images was performed using 
SExtractor (Bertin & Arnouts 1996). SExtractor was 
run in dual image mode using the narrow band image as the 
detection image for each of the broad band images and the 
narrow band image itself. Photometric zero-points were cal¬ 
culated using the narrow-band observations of LTT 9491. As 
there is no NB497 narrow band standard photometry avail¬ 
able for the standard star, we assumed the narrow band 
magnitude to be equal to the quoted V band magnitude of 
V = 14.06 (the B — V colour of this star is +0.03). The zero- 
point magnitude for our observations given a 1 s exposure 
was Mzp = 24.34, which equates to Mzp = 32.04 for our 
1200 s exposures. Taking into account atmospheric absorp¬ 


tion for each field, the airmass ranged between 1.1 — 1.5. 
Taking an extinction coefficient of k = 0.12 gives a range of 
extinction corrections of ~ 0.13 — 0.19. The resulting indi¬ 
vidual zero-points are given in Tab. 1. 

Aperture magnitudes were measured using an aperture 
diameter of 3.0'' (approximately twice the seeing EWHM 
of the broad band data). Total magnitudes were measured 
using MAG_AUTO. la magnitude depths were estimated 
using the errors calculated in 3.0" diameter apertures in 
SExtractor. 

These initial catalogues were used for the selection of 
targets for the VLT VIMOS observations described in the 
following section. 


2.3 Improved Astrometry, Object Detection and 
Photometry 

The initial astrometric solution was discovered to be inac¬ 
curate during the analysis of the VIMOS spectroscopic data 
(as described in Sec. 3). We therefore re-calibrated the as¬ 
trometry and produced improved photometric catalogues for 
all of our fields. 

This re-calibration of the astrometric solution was per¬ 
formed using SCAMP (Bertin 2006) with each stacked narrow 
band image in conjunction with the associated VLRS broad 
band images, thus providing a consistent astrometric solu¬ 
tion across the narrow and broad band images. Since we 
applied the geometric transformation to all images, the ob¬ 
ject positions in each broad-band image were matched with 
objects in the NB band. The distortion correction in the re¬ 
duction process corrected the geometric distortion and we 
obtained good astrometry to better than ±0.2" rms over 
the image. Astrometry was made based on USNO at ESO 
catalogue and ~ 1000 stars identified in the stacked images. 
The position of USNO objects were approximately uniformly 
distributed over the entire stacked images. The absolute co¬ 
ordinates of our objects were obtained from these USNO 
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objects. The images were then transformed based on the 
SCAMP solution using swarp (Bertin et al. 2002). 

We then produced two sets of catalogues based on these 
distortion-corrected images, both sets using SExtractor in 
dual image with one set using the narrow-band images as 
the detection image and the second set using a combined 
image as the detection image. The combined image was 
produced using swarp and in each field combines the narrow 
band and B and R (or V) images in an optimal way for 
source detection (see Bertin et al. 2002, for details). 


3 SPECTROSCOPIC OBSERVATIONS 

3.1 Overview 

We have made spectroscopic follow-up observations of the 
LAE candidates with the Visible Multi-Object Spectrograph 
(VIMOS) (Le Eevre et al. 2003) on the Very Large Tele¬ 
scope (VLT). The VIMOS focal plane is divided into 4 quad¬ 
rants, each measuring 7' x 8^ (with gaps of ~ 2' between 
adjacent detectors). We observed with the HR_Blue grism, 
allowing us to observe ~ 20 — 30 objects per detector with 
a spectral coverage of 4150 < A < 6000 and a resolution of 
R ~ 2,050 (given our 1" slit-width). 

VIMOS Mask Preparation was made by using VIMOS 
Mask Preparation Software (VMMPS)^ and the observa¬ 
tions were conducted at the end of August to the begin¬ 
ning of September, 2011 (part of the observing run ESO-ID 
086.A-0520 B, P.I. H. Erancke). 

The observations were taken in 3 of our 5 LAE fields: 
QSO B2359+068; QSO J0124+0044; and LBQS 0301-0035. 
Details of the observations are given in Table 3. 


3.2 Initial photometric target selection 

LAEs were selected for the VLT VIMOS spectroscopic ob¬ 
servations using the narrow band 497 nm and broad band 
B, V and R band data already described. The transmission 
curves for each of these are shown in Eig. 1. Eor this selec¬ 
tion we used a narrow band magnitude cut of mNB 497 ^ 26 
in combination with (B — mNB 497 ) and {V — ?tinb 497 ) or 
{R — mNB 497 ) colour cuts. These colour cuts were tailored 
individually to each field in order to optimise the slit alloca¬ 
tions in the VIMOS fields targeted. The colour cuts used in 
each field are given in Table 3, along with the total number 
of objects selected (across the whole ~ 0.5° x 0.5° field in 
each case) and the number of these that were targeted in 
the VIMOS masks (i.e. over an area of ~ 16^ x 16^). 

Different groups use different selection criteria to 
search for LAE candidates photometrically. Eor example: 
narrow-band and combined broad-band observations (e.g. 
(B + R)/2); one narrow-band filter and one broad-band 
filter; and one narrow-band filter and multiple individual 
broad-band filters. We choose the latter approach, using a 
two-colour approach which applied R — NB497 colour and 
B — NB497 colour, similar to those in Eynbo et al. (2003); 
Nilsson et al. (2007); Ouchi et al. (2008). 

Eollowing Ouchi et al. (2008), we started from a colour 
cut of (V—mNB 497 ) = 1.2, which equates to a cut in constant 
rest-frame equivalent width of EWo ~ 45 A (or EWo « 65 A 
when taking the Lya Gunn-Peterson trough into account - 
see Ouchi et al. 2008). Eor our fields containing R band data, 
but no V band data, we took into account the mean {V — R) 
colour difference in the 22 < mNB 497 < 25 galaxy population 
giving an equivalent R-band cut is (R — mNB 497 ) = 0.8. In 
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Figure 1. The transmission curve of the filters. The solid line 
indicates the narrow-band filter, NB497. The dashed lines indicate 
broad-band filters, B, V, and R from left to right respectively. 


addition to this, we also applied a cut using the B band of 
B —NB497 > 1.3. In order to maximise the slit allocations in 
the VIMOS masks, we then allowed these colour constraints 
to be relaxed. The resulting colour cuts used for the VIMOS 
target selection in each field are given in Tab. 3. These colour 
cuts were combined with a magnitude limit of mNB 497 ^ 26. 


3.3 Data reduction 

The reduction of the spectroscopic data was performed with 
the VIMOS pipeline esorex packages^. The main procedure 
includes creating master calibration data, reducing science 
frames, and extracting objects. Eollowing the pipeline man¬ 
ual, we firstly created a master-bias with the recipe vmbias. 
An output master-bias, was then used in the reduction of 
the flat field, arc lamp, and scientific exposures. The next 
step was using the recipe vmmoscalib to process flat field ex¬ 
posures and arc lamp exposures. We used the vmmosscience 
recipe to process science frames with the cosmic ray rejec¬ 
tion applied at this point in the process. We then combined 
the images for each field using the iraf imcombine task. 

The object extraction is made by applying an optimal 
extraction algorithm (Horne 1986). The wavelength calibra¬ 
tion was performed using the input wavelength calibration 
and sky lines, and sky background subtracted. We found di¬ 
agonal stray light interfered with the wavelength calibration 
and sky subtraction in quadrants 3 and 4. To fix this, we 
isolated the brightest part of this light from the flat lamp, 
fitted the smooth pattern of the flat and subtracted it out. 


3.4 LAE identification 

We found significant emission within the wavelength range 
of the narrow band filter in 8 of the 23 candidates in the QSO 
B2359-1-0653 field, 13 of the 45 in the QSO J0124-h0044 field, 
and 14 of the 46 in the LBQS 0301—0035 field. These are 
listed in full in Tab. 4, whilst the reduced ID spectra (with 
2D spectra inset) are shown in Eig. Al, Eig. A2 and Eig. A3 
for the three separate fields. In each figure, the solid black 
line shows the reduced ID spectrum, the shaded grey region 
shows the normalised NB497 filter response and the inset 
image shows the 2D spectrum within ±60 A of the peak of 
the NB497 filter response curve. 


^ http://www.eso.org/sci/observing/phase2/SMGuidelines/VMMPS.hHnhttp://www.eso.org/sci/software/cpl/esorex.html 
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Table 3. Details of the VIMOS LAE observations. 


Field 

R.A. Dec. 

(J2000) 

Constraints 

.^LAE 

(Sel/Obs) 

Seeing 

(") 

Integration time 
(hours) 

QSO B2359+068 

00:02:11.457 

+07:15:32.25 

{R — mNB497) > 0-2 
(B — mNB497) > 0-9 

304/37 

1.0 

2.0 

QSO J0124+0044 

01:24:36.236 

+00:51:07.19 

{V — mNB497) > 1-1 
{B — mNB497) > 1-4 

141/45 

1.0 

4.3 

LBQS 0301-0035 

03:03:10.208 

-00:16:20.98 

{R — mNB497) > 0-8 
{B — mNB497) > 1-4 

214/46 

1.0 

3.3 


Table 4. Galaxies with confirmed Lya emission from the VLT VIMOS observations. The top, middle, and bottom entries correspond 
to the galaxies in the fields of the quasars B2359+0653, J01214+0044, and Q0301—0035 respectively 


ID 


RA DEC 

(J2000) 

z 

771497 

B - m497 

(AB) 

R — 777497 

VLRS J000223.50+071321.0 

0.597932 

7.222500 

3.0939 

23.48 

2.44 

1.50 

VLRS J000243.85+071021.5 

0.682742 

7.172648 

3.1079 

24.49 

2.75 

4.40 

VLRS J000222.05+071010.6 

0.591878 

7.169614 

3.0543 

24.81 

1.96 

2.91 

VLRS J000221.03+070915.6 

0.587649 

7.154351 

3.1187 

24.69 

2.11 

3.34 

VLRS J000150.11+072144.3 

0.458817 

7.362311 

3.1112 

25.57 

— 

— 

VLRS J000200.93+071939.3 

0.503911 

7.327593 

3.0994 

25.46 

— 

1.43 

VLRS J000143.31+071209.4 

0.430497 

7.202622 

3.1015 

25.55 

— 

1.40 

VLRS J000154.47+070912.2 

0.476987 

7.153391 

3.1052 

24.17 

3.99 

2.28 

ID 


RA 

DEC 

z 

777497 

B - 7T7497 

V - 777497 



(J2000) 



(AB) 


VLRS 

J012457.44+004648.0 

21.239342 

0.780015 

3.1030 

24.70 

1.68 

3.29 

VLRS 

J012451.71+004502.5 

21.215461 

0.750693 

3.0520 

25.18 

2.15 

1.30 

VLRS 

J012450.74+004352.3 

21.211416 

0.731202 

3.0811 

24.55 

— 

0.24 

VLRS 

J012445.90+004317.0 

21.191265 

0.721403 

3.0706 

25.23 

1.81 

1.12 

VLRS 

J012509.24+005710.7 

21.288526 

0.952986 

3.1228 

24.84 

3.05 

— 

VLRS 

J012503.13+005634.7 

21.263065 

0.942975 

3.1104 

23.71 

3.22 

3.21 

VLRS 

J012408.25+005656.9 

21.034377 

0.949161 

3.0748 

24.15 

3.02 

4.67 

VLRS 

J012416.71+005622.2 

21.069636 

0.939499 

3.0672 

25.20 

3.21 

4.52 

VLRS 

J012413.08+005529.5 

21.054527 

0.924861 

3.0734 

25.54 

0.44 

0.86 

VLRS 

J012408.69+005341.3 

21.036244 

0.894817 

3.0930 

24.94 

2.60 

3.74 

VLRS 

J012428.28+005249.8 

21.117851 

0.880499 

3.0761 

24.00 

— 

1.53 

VLRS 

J012406.88+004750.3 

21.028689 

0.797311 

3.0792 

24.39 

2.04 

2.39 

VLRS 

J012417.96+004733.7 

21.074841 

0.792696 

3.0786 

24.25 

3.22 

3.07 

ID 


RA 

DEC 

z 

777497 

B — 777497 

R - 777497 



(J2000) 



(AB) 


VLRS 

J030326.84-001805.8 

45.861860 

-0.301630 

3.0951 

25.96 

— 

0.14 

VLRS 

J030337.40-001959.3 

45.905844 

-0.333158 

3.1322 

25.67 

0.54 

0.39 

VLRS 

J030319.11-002139.7 

45.829661 

-0.361044 

3.0646 

25.89 

— 

— 

VLRS 

J030325.62-002350.5 

45.856789 

-0.397361 

3.0757 

25.28 

2.96 

— 

VLRS 

J030326.31-002400.2 

45.859633 

-0.400077 

3.0541 

24.82 

1.48 

1.05 

VLRS 

J030333.31-001047.6 

45.888794 

-0.179914 

3.0746 

24.57 

2.07 

— 

VLRS 

J030256.04-001007.2 

45.733499 

-0.168692 

3.0802 

24.32 

2.47 

2.51 

VLRS 

J030248.36-001022.7 

45.701519 

-0.172993 

3.0943 

25.52 

— 

-0.01 

VLRS 

J030241.28-001140.6 

45.672003 

-0.194622 

3.1031 

25.28 

3.09 

— 

VLRS 

J030244.58-001341.8 

45.685771 

-0.228293 

3.0730 

23.87 

1.52 

1.21 

VLRS 

J030254.08-001414.6 

45.725362 

-0.237397 

3.0854 

24.08 

3.67 

2.24 

VLRS 

J030244.33-001911.9 

45.684716 

-0.319997 

3.1251 

24.37 

1.41 

0.62 

VLRS 

J030257.61-002058.5 

45.740049 

-0.349601 

3.0939 

24.76 

— 

3.26 

VLRS 

J030257.31-002301.0 

45.738794 

-0.383629 

3.1019 

25.62 

1.90 

0.82 


Given all our detections are single line detections, the 
possibility exists that some of these could be [Oil] emission 
from faint z ~ 0.33 galaxies. The spectral resolution of our 
observations is marginally sufficient to discern the double 
peaked emission inherent in the Oil double lines (the reso¬ 
lution is ~ 2.5 A while the [Oil] doublet is about 2.7 A). 
Although the resolution is very close to the peak separa¬ 
tion, even in the event of the doublet being smoothed out 
in the spectrum by the instrument response, the presence 
of the [Oil] doublet would lead to a significantly broad line. 
We have analysed each emission line, plotting the doublet 


separation over each and evaluating the likelihood of each 
line being [Oil] emission. In most cases (> 90%), clearly 
shows no sign of being a doublet (i.e. there is no double 
peak and the single peak emission is not broad enough to 
be the two [Oil] lines convolved into one via the instrument 
response). For the remaining < 5%, we cannot discount the 
possibility of the emission being low redshift [Oil], but at 
the same time, none are clear-cut cases of [Oil] emission. 
These ambiguous < 5% are broad emission that show no 
significant detections of double peaks within the spectral 
noise. Indeed higher signal-to-noise and resolution observa- 
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tions would be needed to fully discern any double peaked 
nature to these broader lines. In addition, we note that even 
if double-peaked emission were detected, the complex nature 
of the escape of Lya photons from galaxies can also produce 
double emission peaks. In all we find no strong evidence of 
any double emission in any of the detected emission lines 
and conclude that > 95% of these lines are Lya emission at 
z 3.06. 

The success rate of the observations were somewhat low 
in part due to the relaxed constraints used to select candi¬ 
dates. Additionally, the slit mask for the second quadrant in 
the QSO B2359+0653 field was not properly aligned in its 
mount leading to no detections in that quadrant. However, 
in addition to this, a poor astrometric solution was found 
to have contributed to the low success rate. The astromet¬ 
ric errors for the images created two issues: 1. in the QSO 
B2359+0653 field, a number of targets were not correctly 
aligned within the slits during the observations; and 2. in¬ 
consistencies between the astrometry between the narrow 
band and broad bands led to incorrect broad band photom¬ 
etry measurements. Apertures placed to measure the pho¬ 
tometry in the broad band images were thus offset from the 
intended targets in some cases, leading to incorrectly faint 
magnitudes being measured and objects with no emission 
being promoted into the selection criteria. 

Following the spectroscopic observations, we thus re¬ 
calculated matched astrometric solutions for all the data as 
described in Sec. 2.3. 


3.5 Selection efficiency and optimised selection 
criteria 

The selection success rates are relatively low in our spectro¬ 
scopic sample of LAE candidates (31%). This is the result 
of the original criteria being comparatively flexible in order 
to allow the maximal number of slits placed per VIMOS 
field, whilst as discussed the astrometric calibrations were 
sub-optimal prior to the spectroscopic follow-up, causing er¬ 
rors in the photometry. In this section we therefore provide 
an optimised and uniform set of selection constraints, max¬ 
imising the numbers of successfully identified LAEs based on 
our updated photometric catalogues and the spectroscopic 
observations. 

First of all, we now use a narrow band magnitude cut of 
TnNB 497 ^ 25.0 (as opposed to mNB 497 ^ 26 previously), as 
in all three fields there is significant scatter from the galaxy 
locus at fainter magnitudes than this due to the measure¬ 
ment uncertainties on the photometry. The converse of this 
is that the measured colours of many LAEs beyond this limit 
are scattered into the galaxy locus in the colour-magnitude 
diagram. 

Secondly, we optimise the E —mNB 497 , —mNB 497 and 

B — mNB 497 colour cuts at ?tinb 497 ^ 25.0 in order to max¬ 
imise the numbers of LAEs found and minimise the num¬ 
ber of interlopers given the spectroscopic observations. We 
also apply the constraint that the three colour criteria cor¬ 
respond to the same equivalent width constraint. For this 
we simulate LAE spectra by assuming a /3 = —1.4 power- 
law slope (e.g. Berry et al. 2012) combined with a range of 
equivalent Lya emission lines. We find that taking colour 
constraints of E — mNB 497 ^ 1.2, R — mNB 497 ^ 1.0 and 
B — mNB 497 ^1.4 selects 18 out of the 20 spectroscopically 
confirmed ?tinb 497 LAEs at mNB 497 ^ 25.0, whilst missing 2 
of these. The selection also includes 5 non-detections from 
the spectroscopic observations that may be contamination or 
weakly emitting z ~ 3.1 LAEs (likely scattered into the se¬ 
lection region due to noise on the photometry). These colour 
cuts are equivalent to a rest frame equivalent width cut of 


Table 5. Number of LAE Candidates from our selection in each 
field. 


Field 

^LAE 

Plae 

(arcmin“^) 

QSO B2359+068 

124 

0.160 

QSO J0124+0044 

139 

0.186 

LBQS 0301-0035 

170 

0.201 

PKS 2126-158 

140 

0.158 

QSO B2231-0015 

70 

0.119 

All fields 

643 

0.167 


~ 45 A, when assuming the simple /3 = —1.4 continuum 
power law slope. We do not however include the effect of 
the Gunn-Peterson trough and so, as discussed and mod¬ 
elled in Ouchi et al. (2008), our nominal 45 A selection 
corresponds to a 65 A cut when taking this into account. 

Additionally, we also reject objects with S Extractor 
flags greater than zero; mask low signal-to-noise regions of 
the images (primarily the image edges) to avoid noise con¬ 
tamination in these regions; and (in order to minimise the 
contribution of artefacts in the images) apply elongation 
(esEx) and FWHM limits. By inspection of the images, ob¬ 
jects with a measured elongation of esEx > 2.8 are exclu¬ 
sively artefacts (primarily cosmic ray hits not removed by 
crreject). In addition, objects with measured EWHM less 
than the measured image FWHM are also rejected (on in¬ 
spection these are also exclusively artefacts and primarily 
remnant cosmic rays). 

The selection criteria are summarised as follows: 

• 20 < mNB497 < 25; 

• (E - mNB497) > 1.2 or (R - mNB497) > 1.0; 

• {B — mNB497) > 1.4; 

• FLAGsex = 0; 

• esEx < 2.7; 

• FWHM ^ FWHMi^age. 

These colour cuts are shown by the solid lines in Figs. 2 
and 3. In each figure the galaxy population is shown by the 
2d histogram, whilst selected LAE candidates are shown by 
open circles. The dashed horizontal lines show the position 
of the locus in (E —mNB 497 ) or (i^ —mNB 497 ) colour (Fig. 2) 
and {B — rriNBAor) colour (Fig. 3). The grey and black filled 
circles show LAE candidates observed with VLT VIMOS 
and are described further in Sec. 3. The seeing for each 
field (i.e. FWHMimage) is given in Tab. 1. As discussed, we 
have produced catalogues based on object detection in the 
narrow-band images alone and also in a combined image 
of the narrow and broad band images. We have run our se¬ 
lection on both of these sets of catalogues. In the results that 
we report here, we use the y^-image detected catalogues. We 
have run all our analysis on the narrow-band image detected 
catalogues and find consistent results within the estimated 
uncertainties. 

The numbers and sky densities of candidates in each 
field given by these criteria are given in Tab. 5. Taking only 
those spectroscopically observed targets that fall within the 
updated selection criteria, we find a success rate of 78±18%. 
The proportion of targets selected in this new sample that 
were also present in the original target selection is ~ 40%, 
thus calculating the success rate for our updated selection 
using the spectroscopic sample is arguably only represen¬ 
tative of 40% of the sample. However, the remaining 60% 
of targets selected in the updated selection were missed in 
the original selection due to astrometric errors. The loss of 
targets due to the sub-optimal astrometric solution has no 
connection with intrinsic source properties (rather it is a 
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Figure 2. Colour magnitude diagrams showing V — mNB497 versus mNB497 for the QSO J0124+0044 field and R — mNB497 versus 
mNB 497 for the LBQS 0301-0035, PKS 2126-158, QSO B2231-0015, and QSO B2359+068 fields. The shaded 2D histogram shows 
the general galaxy population. Grey and black filled circles show LAE candidates observed with VLT VIMOS, the grey showing those 
with no spectroscopic detection and the black showing those confirmed as 2 ; = 3.06 LAEs. The region to the upper left of the solid 
horizontal and vertical bounding lines give the updated selection region, whilst the open circles show LAE candidates selected based on 
these bounds. Note that the objects identified as being targeted with VIMOS that lie within the main galaxy population are those that 
were subject to significant photometric errors due to the astrometric mis-alignment of images (see Sec. 3.4). 


function of the on-sky coordinates of the sources) and so 
should not introduce any biases between the 40% that were 
selected and the 60% that were not, so it would be reasonable 
to assume that the result for the 40% that were included, 
holds for the 60% that were missed due to astrometric issues. 

Comparing to the success rate in other’s work, 
Ouchi et al. (2008) identified line emitters from 60% of their 
targets at this redshift, whilst Fynbo et al. (2001, 2003) re¬ 
ported the spectroscopic follow-up success rate of 75 - 90% 
for z ~ 3 LAE surveys. 

The distribution of the selected LAE candidates in each 
field is shown in Eig. 4, with the same symbols as in previous 
plots (i.e. open circles for LAE candidates and filled grey and 
black circles for candidates observed with VLT VIMOS). 
Eilled stars in Eig. 4 denote the positions of background 
QSOs in the fields. 

3.6 Redshift and spatial distribution of LAEs 

The redshift distribution, n(z), for the VLT VIMOS con¬ 
firmed LAEs is shown in Eig. 5 (solid blue histogram). Eor 
comparison, the throughput curve of the NB497 filter is also 
shown (solid black curve). 

The survey sky coverage is illustrated in Eig. 4, where 
the distributions of candidate LAEs (small grey circles) and 
confirmed LAEs (large black circles) are shown for each 
field. Background z > 3.1 QSOs are also shown (filled- 
stars). These five fields cover a total area (after masking) 
of 1.07 deg^. 


Integrating along the normalised VIMOS redshift dis¬ 
tribution and combining it with this total sky area gives 
an LAE survey volume of 2.51 x 10^ /i“^Mpc^. Integrating 
instead the NB497 filter response curve gives a somewhat 
higher survey volume estimate of 3.35 x 10^ h~^Mpc^. In 
the analysis that follows, we take the volume given by the 
filter response curve given that the measured redshift distri¬ 
bution is somewhat under-sampled. Eor reference, a top-hat 
redshift distribution (between z = 3.05 and z = 3.14) gives 
a survey volume of 4.69 x 10^ h~^Mpc^. 

3.7 Selection Completeness 

In order to estimate the selection completeness, we first es¬ 
timate the percentage completeness as a function of source 
magnitude of the individual survey images. To do so, we 
simulate sources in the images in set magnitude bins. We 
attempt to detect the simulated sources using SExtractor, 
identically to the extraction process for the original data. Eor 
each magnitude bin, we simulate 100 sources, with the num¬ 
ber recovered giving the fractional completeness in each bin. 
We avoid source confusion by only placing simulated sources 
in empty regions with no genuine sources found within 2.5x 
the image EWHM. 

We applied these detection rates with simulated LAE 
spectra with a range of equivalent widths and assuming a 
spectral slope of (3 = —1.4 (e.g. Berry et al. 2012). The re¬ 
sults of this simulation for each field are shown in Eig. 6. 
Each panel shows one of our five individual fields, with the 
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Figure 3. Colour magnitude diagrams showing B —mNB497 versus mNB497 for our five QSO fields: QSO J0124+0044, LBQS 0301—0035, 
PKS 2126—158, QSO B2231—0015, and QSO B2359+068. The shaded 2D histogram shows the general galaxy population. Grey and 
black filled circles show LAE candidates observed with VLT VIMOS, the grey showing those with no spectroscopic detection and the 
black showing those confirmed as 2 : = 3.06 LAEs. The region to the upper left of the solid horizontal and vertical bounding lines give the 
updated selection region, whilst the open circles show LAE candidates selected based on these bounds. Note that the objects identified 
as being targeted with VIMOS that lie within the main galaxy population are those that were subject to significant photometric errors 
due to the astrometric mis-alignment of images (see Sec. 3.4). 


shaded regions showing the range of our LAE selection in 
R-EW space. The contour curves show the 20%, 40%, 60%, 
and 80% completeness levels. 

Overall, we find good consistency across our five fields 
in terms of the completeness profile. At fainter continuum 
luminosities (i.e. R > 25.5), we find that we increasingly 
miss the high equivalent width range of our selection with 
all 5 fields falling to 20% completeness at R ^ 27. 

4 NUMBER COUNTS LUMINOSITY 

FUNCTION 

4.1 Narrow band number counts 

Using the updated selection criteria given in the previous 
section, we calculate the galaxy number counts for the LAE 
sample as a function of narrow band magnitude. This is 
shown in Eig. 7, where the filled black circles show our re¬ 
sults. The error bars show field-to-field error estimates on 
the points. Eor reference, we also show the NB497 number 
counts for all sources in our fields (grey circles). 

Eor comparison, we show the z ~ 3.1 narrow band num¬ 
ber counts published by Gronwall et al. (2007, diamonds) 
and Ouchi et al. (2008, triangles). Gronwall et al. (2007) 
conducted a 0.28 deg^ survey at 2 : ~ 3.1 in the Extended 
Ghandra Deep Eield South (EGDES). Their selection is 
based on a measured colours derived from the NB5000 fil¬ 
ter and a broad band image constructed using a combi¬ 
nation of images taken using the B and V band filters. 


Their criteria were thusimNBsooo — tub+v < 1.03, with a 
narrow-band magnitude limit of ttinbsooo ^ 25.4. This gave 
a sample of 162 galaxies with rest-frame equivalent widths 
of > 20 A (approximately 80 A in the observer’s frame). The 
Ouchi et al. (2008) selection criteria are based on a cuts of 
U —ttinbsos >1.2 and B — V > 0.5. They detected 356 pho¬ 
tometric LAEs, with rest-frame equivalent widths of ^ 65 A, 
over a ~ 1 deg^ area. 

A direct comparison of the Gronwall et al. (2007) and 
(Ouchi et al. 2008) number counts of 2 : ~ 3.1 LAEs is 
shown by Giardullo et al. (2012). Giardullo et al. (2012) ex¬ 
tended the study of Gronwall et al. (2007) by re-imaging the 
EGDE-S with a 57 A EWHM nearly top-hat filter centred 
at 5010 A. They found a total number of 360 z ^ 3.1 
LAEs, a subset of which are also selected by Gronwall et al. 
(2007). They then reproduced the luminosity function and 
concluded that their result is statistically identical to values 
in Gronwall et al. (2007) and Ouchi et al. (2008), once dif¬ 
ferences in the filter EWHM and EW selection constraints 
were accounted for. 

The VLRS LAE number counts are consistent with the 
number counts of Gronwall et al. (2007) at the ~ la level, 
and marginally higher than the counts of Ouchi et al. (2008) 
(although still only at the ^ 1 — 2a level). Our plotted 
number counts include a factor of 0.78 to account for the 
number of non-detections when our sample was observed 
spectroscopically. This is perhaps a conservative estimate 
for the success of our selection criteria in identifying LAEs, 
given that some fraction of those non-detections could have 
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Figure 4. The distribution of candidate (grey filled circles) and confirmed (black filled circles) LAEs in the five survey fields. Known 
QSOs at the LAE redshifts or higher (i.e. 2 ; > 3.05) are also plotted for reference (stars). The dashed green boxes show the areal coverage 
of the VIMOS quadrants for the spectroscopic observations (note that the fields PKS 2126-158 and QSO B2231-0015 fields have not been 
observed spectroscopically). The central QSO in the QSO B2231-0015 field is not plotted as it has a redshift lower than that of the LAE 
selection (2 ;qso = 3.027). 



Figure 5. Redshift distribution for LAE spectroscopic identifi¬ 
cation from VIMOS (normalised histogram) and the NB497 filter 
with wavelength converted to redshift (solid curve). The NB497 
filter is normalised such that the peak is unity. 

been LAEs below the detection threshold of the spectro¬ 
scopic observations. In the most extreme case, i.e. if all the 
non-detections were in fact genuine LAEs, the VLRS LAE 
number counts shown in Fig. 7 could be increased by a fac¬ 
tor of 25%. This would perhaps cause some tension with 
the counts of Ouchi et al. (2008) at the 2 — 3cr level. 

For context in comparing the individual results, the 
equivalent width limits for each survey are: 65 A (this pa¬ 
per), 20 A (Gronwall et al. 2007), and 65 A (Ouchi et al. 


2008). Ciardullo et al. (2012) suggest that the lower equiv¬ 
alent width limit of Gronwall et al. (2007) should measure 
10% higher numbers of LAEs compared to Ouchi et al. 
(2008), and hence our own sample. In addition to the 
different EW cuts, given the narrower filter used by 
Gronwall et al. (2007), one would expect to measure lower 
counts by a factor of | given the smaller volume coverage 
this entails. Giardullo et al. (2012) do not explicitly state 
that the difference in filter width is included in their predic¬ 
tion of a difference of 10% in the numbers of LAEs between 
the samples, however their estimate would seem somewhat 
low if it is not (for example when compared to the equiva¬ 
lent width distribution of LBGs measured by Shapley et al. 
2003). 

Although our main aim has been to check the represen- 
tiveness of our LAE counts as a prelude to using the LAE 
for clustering studies, we now briefly compare our counts 
to a recent theoretical model. Shown in Fig. 7 therefore 
is the Orsi et al. (2012) model for LAEs updated to the 
latest variant of the GALFORM semi-analytical model of 
galaxy formation (Lacey et al. 2015). This model for LAEs 
makes use of a Monte Garlo Lya radiative transfer code to 
compute the escape fraction of LAEs at different redshifts. 
The GALFORM model is calibrated to match a number of 
observational data, mostly at z = 0 (e.g. optical and NIR 
LFs, Mbh — Mbuige relation, morphological fractions), but 
also sub-mm, FIR counts and UV luminosity functions at 
high redshift. A radiation transfer model was then used to 
model the Lya escape fraction, based on an outflow model 
which is itself dependent on the individual galaxy proper¬ 
ties. This radiation transfer and outflow model (which in 
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Figure 6. Estimated selection completeness in the five observed 
fields based on simulated sources in the survey images. The 
shaded region shows the extent of our selection criteria in i?—EW 
space, whilst the contour curves show the 20%, 40%, 60%, and 
80% completeness levels. 



Figure 7. The number counts of photometrically detected LAEs 
based on the selection criteria given in Sec. 3.5 (filled black cir¬ 
cles). The diamonds and triangles show number counts for com¬ 
parable LAE selections by Gronwall et al. (2007) and Ouchi et al. 
(2008) respectively. Small filled grey circles show the number 
counts of all objects detected in the narrow band imaging. The 
solid curve shows model LAE number counts from the GAL- 
FORM galaxy formation model. The shaded curve shows the LAE 
R-band counts adapted from Fig. 8 with a shift of Am = —1.6 to 
overlay the narrow band counts. 
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Figure 8. The R-band luminosity function of the VLRS 2 ; = 3.1 
LAEs (filled circles). For comparison, we also show the equiv¬ 
alent LAE luminosity functions from Gronwall et al. (2007, di¬ 
amonds) and Ouchi et al. (2008, triangles). Also shown is the 
observed r-band luminosity functions of 2 ; ~ 3 LBGs from 
Reddy & Steidel (2009, pentagons) and from the GFHT fields of 
van der Burg et al. (2010, stars). The dotted and dash-dot curves 
show the Reddy & Steidel (2009) luminosity function reduced by 
a factor of 3 and a factor of 10 respectively, whilst the dashed 
curve shows the GALFORM model result. 


this case consists of an expanding thin shell of material) 
was tuned to match the luminosity function of Ouchi et al. 
(2008). Given the agreement between our own results and 
those of Ouchi et al. (2008), it is no surprise to find the 
model successfully reproduces our own narrow band number 
counts at mNBA 97 ^ 24.2, although the agreement worsens 
at fainter magnitudes (as it does with the Ouchi et al. 2008 
results). Even to the extent that they agree, the tuning of 
the model to the luminosity function of Ouchi et al. (2008) 
means that the observed LAE number counts are therefore 
not a test of the model. This and the disagreement at faint 
fluxes must be borne in mind when comparing the LAE clus¬ 
tering data with the model in Sec. 5.2. 


4.2 Continuum luminosity function 

We estimated the i7-band continuum luminosity function 
of our sample of LAEs. As in Gronwall et al. (2007), we 
use the filter curve to define the survey volume used for 
the luminosity function calculation, which gives a volume 
of 3.35 X 10^ h ^Mpc^. The resulting continuum luminosity 
function evaluated across all 5 fields (and scaled in accor¬ 
dance with the estimated 78% success rate of our sample) is 
shown in Eig. 8 (filled black circles). Plotted points are cor¬ 
rected for the incompleteness as a function of R magnitude 
estimated in Sec. 3.7. The error bars give the uncertainty 
based on field-to-field estimates using our five imaging fields 
in combination with the uncertainty introduced by the esti¬ 
mate of the same success rate. 

We show our results in comparison to those of 
Gronwall et al. (2007, diamonds) and Ouchi et al. (2008, tri¬ 
angles) . Our results are consistent with those of Ouchi et al. 
(2008) within the calculated errors to a magnitude of 
R ^ 26. In contrast, our luminosity function (and that 
of Ouchi et al. 2008) is only consistent with that of 
Gronwall et al. (2007) at magnitudes of > 25. At < 25, 
the Gronwall et al. (2007) luminosity function is signifi¬ 
cantly elevated above our own and approaches the vol¬ 
ume densities of LBGs at the same redshift - shown by 
the filled stars (van der Burg et al. 2010) and pentagons 
(Reddy Steidel 2009). Indeed, whilst Ouchi et al. (2008) 
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report their LAE luminosity function to be ~ i of that of 
the Reddy & Steidel (2009) LEG luminosity function (rep¬ 
resented by the dash-dot curve in Fig. 8), Gronwall et al. 
(2007) find that their LAE Muv luminosity function is 
equivalent to I of the Reddy & Steidel (2009) LEG lumi¬ 
nosity function (represented by the dotted curve). Our own 
LAE Muv luminosity function is consistent with LAEs (of 
equivalent width > 65 A) being lOx less common than 
LEGs, given magnitudes of R < 26 (Muv ^ —19.6). The 
difference between the Gronwall et al. (2007) and our own 
results would appear to be driven by the differing equiv¬ 
alent width cuts used. This is supported by the distribu¬ 
tion function of equivalent widths in LEGs as presented 
by Shapley et al. (2003), whereby there are ~ 3x as many 
galaxies with equivalent widths in the range > 20 A (equiv¬ 
alent to the Gronwall et al. 2007), than in our range > 65 A 
at magnitudes of R < 25.5. 

We note that as with the number counts in the previous 
section, we have applied a correction on the assumption that 
our sample successfully identifies LAEs in 78% of cases. In 
the case of this being overly pessimistic, the presented lumi¬ 
nosity function could be boosted by up to ~ 25%. Within 
the plotted error bars, this does not change any of the com¬ 
parisons discussed above. 

Again, we also present the GALEORM model results 
(which are calculated based on our own selection criteria), 
shown by the dashed curve in Eig. 8. Our observations and 
the model are consistent within the observational uncertain¬ 
ties at R < 25, however the model under-predicts the num¬ 
bers of galaxies at magnitudes fainter than this limit (con¬ 
trasting with the narrow band counts in which the model 
over-predicts the faint end). We note that, when the specific 
selection constraints of Gronwall et al. (2007) are applied to 
the GALEORM model, the model under-predicts the contin¬ 
uum luminosity function at both the bright and faint ends, 
whilst when the same is done with the Ouchi et al. (2008) 
selection constraints, the result is the same as seen with our 
own data. These comparisons between our results and pre¬ 
vious observational results and the model will be instructive 
when analysing the different clustering results. 

Returning to Eig. 7, we show the continuum number 
counts (shifted by Am = —1.6, pale red shaded curve) in 
comparison to the narrow band number counts. In both the 
narrow band and continuum number counts (and the con¬ 
tinuum luminosity function), we see evidence for a turn over 
at faint magnitudes. Indeed, Ouchi et al. (2008) for exam¬ 
ple measure ~ 5 x 10^^ ergs s“^ for their z 3.1 

LAEs, which corresponds to a narrow band magnitude of 
mNB ~ 24.5. Our results support the assertion that such 
LAE samples reaching depths of mNB ~ 25 are probing the 
knee of both the NE and continuum luminosity functions. 
Eriefly comparing the narrow-band and continuum number 
counts, we find that the continuum number counts appear 
steeper, hinting at an evolution in the EW distribution as a 
function of galaxy brightness. 


5 LAE AUTO-CORRELATION FUNCTION 

We now measure the clustering properties of the photometric 
LAE sample to derive halo masses and other properties. 


5.1 Clustering Estimator 

We estimate the angular auto-correlation function, w{0)^ 
using the Landy-Szalay estimator (Landy & Szalay 1993), 
which is given by: 


DD(0)-2DR(0) + RR(0) 

=- wm - ^ ^ 

where DD(^), DR(^), and RR(^) are the numbers of galaxy- 
galaxy pairs, galaxy-random pairs and random-random pairs 
as a function of 0. Eor each field, we generated uniform ran¬ 
dom points with the same area as our masked LAE samples, 
with 20X the number of LAE candidates in each field. Eig. 9 
shows the distribution of the random points used in each 
field for the clustering calculation, showing the extent of the 
fields and the masking implemented. 

We use two error estimators here: simple Poisson er¬ 
ror estimates; and field-to-field error estimates. The Poisson 
estimate is given by: 


crpoi(^) = 


1 + w(0) 

VDD{e)/2 


( 2 ) 


The field-to-field error estimate is given by the error 
on the mean of the measurement across the fields and is 
calculated using: 


crFtF(^) = 


1 


N N 


1 ^ 

— y;[«;i(6») - w(6>)]2 


(3) 


where N is the number of fields (i.e. A = 5), Wi{0) is a 
measurement from the zth field and w{0) is the mean value. 


5.2 Angular Clustering results 

Eig. 10 shows the results of the clustering analysis (filled 
black circles). We use field to field errors everywhere except 
at ^ < 1', where Poisson estimates are used. The plotted 
points include an amplitude correction for contamination in 
the photometrically selected sample, which takes the form: 




Ao 

(1 - fcT 


(4) 


where Ai and Ao are the intrinsic and observed correlation 
amplitude respectively, and fc is the contamination fraction 
(i.e. fc = 0.22). The plotted points have all been adjusted 
upwards by a factor of 1/(1 — /c)^ = 1.64. The points are 
also corrected for the integral constraint, which accounts for 
the effect of finite field sizes. This is estimated following the 
method outlined in Eielby et al. (2014) and takes a value of 

0.024. 

We compare our results in the top panel of Eig. 10 with 
several measurements from other authors. Squares show the 
result from Hayashino et al. (2004) which is obtained from 
283 z = 3.1 LAE candidates with a narrow band magnitude 
limit of mNB 497 <25.8 and EWrest ^ 38 A observed in the 
SSA22a field using the Subaru Telescope. This was a further 
observation of SSA22a of Steidel et al. (2000) who found 72 
LAEs. However, the observed area in Hayashino et al. (2004) 
is about 10 times larger than in Steidel et al. (2000). Eoth 
Steidel et al. (2000) and Hayashino et al. (2004) concluded 
that no significant clustering of z = 3.1 LAEs has been 
shown inside the SSA22a area. We also compared our results 
with Gawiser et al. (2007) (represented by diamonds) who 
studied the clustering properties of 162 LAEs at z = 3.1 at 
a narrow-band completeness limit of msooo ^ 25.4 which 
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Figure 9. The distribution of the random catalogues used in the correlation function analysis, highlighting the masked regions in each 
of the fields. 


were observed in the deep narrow-band MUSYC survey. 
Their EWrest limit is ~ 20 . Ouchi et al. (2010) used 356 
z = 3.1 LAEs from Ouchi et al. (2008) who carried out 
a narrowband survey in the 1 deg^ Subaru/XMM-Newton 
Deep Survery (SXDS) to measure the correlation functionat 
mNB 5003 ^ 25.3 (triangles). 

Comparing our results with each measurement, we 
found that our results agree with Ouchi et al. (2010) within 
their error bars. At larger scales, Hayashino et al. (2004) and 
Gawiser et al. (2007) show slightly higher clustering ampli¬ 
tude than ours, but are also in agreement with our results 
within the error estimates. Given the different equivalent 
width, and flux limits of each survey, this may be evidence 
that the clustering amplitude of LAEs is particularly insensi¬ 
tive to the Lya sample selection. In other words, LAEs trace 
the density field in the same way regardless of their selection. 
Care should be taken when comparing the w{0) measure¬ 
ments of samples with differing filter widths however. Whilst 
Hayashino et al. (2004), Ouchi et al. (2010), and our own 
sample use a filter with EWHM 75 A, the Gawiser et al. 
(2007) LAE sample uses a filter with EWHM 50 AA more 
reliable comparison can be made when comparing the ro 
measurements that follow, which take into account the filter 
profile. Overall, our results provide the most accurate mea¬ 
surement of the LAE clustering signal to date given our large 
total area and the use of multiple fields, which minimise any 
cosmic variance effects. 

In the lower panel of Eig. 10 we compare our results to 
the predictions of the GALEORM model. The filled circles 
are the same as in the top panel, whilst the solid curve and 
shaded regions show the median clustering measured from 
225 mocks and the la and 2cr uncertainties. Each mock cov¬ 
ers a survey area of 1.07 deg^ (i.e. equivalent to our own 
survey area) and so the uncertainties given by the shaded 


regions give a further estimate of the uncertainties in our ob¬ 
servations (and one that corroborates our estimates based 
on the data itself). 

We find that the model results over-predict the ob¬ 
servational measurement at the ~ 2a level. Although the 
semi-analytical model is able to reproduce the narrow-band 
and continuum luminosity function of LAEs at bright mag¬ 
nitudes, it is worth noting that the clustering results will 
be dominated by the more numerous faint LAEs, which 
the model fails to reproduce. It may be that the discrep¬ 
ancy here is connected to the model number count dis¬ 
agreements at faint magnitudes. Taking our results along¬ 
side comparable over-predictions of clustering measurements 
by semi-analytical simulations (e.g. Guo et al. 2011, 2013; 
Bielby et al. 2014), it is clear that clustering measurements 
provide an important additional constraint on galaxy forma¬ 
tion models. Indeed, resolving these discrepancies at faint 
magnitudes in the LAE numbers and clustering remains a 
challenge for the semi-analytical models. 


5.3 Real-space clustering 

We now parameterise the LAE clustering in terms of the 
real-space clustering length as measured via the angular 
two-point correlation function w{0)^ which is effectively a 
weighted projection of the spatial two-point correlation func¬ 
tion ^(r). It is common to transform from the spatial to the 
angular correlation function using Limber’s approximation 
(Limber 1953), however this requires that the depth of the 
galaxy survey, tt (i.e. distance probed along the line of sight), 
is much greater than the on sky maximum separation, a (e.g. 
Limber 1953; Simon 2007). Eor our survey volume, this is not 
the case - the depth of our survey is tt ~ 60 h~^Mpc and 
the width is cr 40 h~^Mpc. We therefore use the full an- 
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Figure 10. Top panel: The angular correlation function of the 
z = 3.1 LAEs. The filled black circles show the result based on our 
VLRS LAE sample. Also shown are model fits to the VLRS LAE 
correlation function based on using the n{z) from our VIMOS 
spectroscopically observed sample (solid curve) and using an n{z) 
based on the NB497 filter profile (dashed curve). A number of 
literature results are also plotted: Ouchi et al. (2010) is shown by 
the blue triangles; the Gawiser et al. (2007) result is shown by 
the orange diamonds; and the result of Hayashino et al. (2004) is 
shown by the cyan squares. Lower panel: The same VLRS LAE 
result (filled black circles) alongside the LAE clustering prediction 
from the GALFORM model. The shaded regions show the 68th 
and 95th percentile uncertainties based on the standard deviation 
of 225 1.07 deg^ mock catalogues. 

alytical form in transforming from the real-space clustering 
form to the projected angular clustering. 

As the underlying model for the real-space clustering, 
we assume a power-law with a slope of 7 = 1.8 (e.g. Peebles 
1980): 

(5) 

where r is the real-space separation between two galaxies 
and ro is the clustering length parameter. 

The spatial correlation function, ^{r) can be related to 
w{0) by (e.g. Phillipps et al. 1978): 

where f(z) is the radial distribution of sources which is given 
by: 

f{z) = X^{z)^^^^nc{z)4>{z) (7) 

where 0 (z) is the selection function of the sample, nc{z) is 
the comoving number density of the sources, x is the ra¬ 
dial comoving distance, and r = r{ 0 ,zi,Z 2 ) is a comoving 
separation between two points at and Z 2 - 


r = Vx^{zi) + X^{z2) - 2 x(zi)x(z2)cos9 (8) 


We assume that over the redshift range probed, the ac¬ 
tual number density of LAEs, nc{z), is constant, whilst we 
have calculated solutions with each of the filter response 
curve and the redshift distribution of the VIMOS detected 
sources as the selection function, 0(^). We performed a min¬ 
imised x^ fit to determine the clustering length, ro, given our 
w{0) clustering measurement and using Eq. 6 to transform 
from the real-space power law correlation function, ^(r), and 
the angular correlation function, w{0). With the n{z) of 
the spectroscopically confirmed LAEs and a fixed slope of 
7 = 1.8, we derived a result of ro = 2.79 zb 0.34 h~^ Mpc. 
The resulting w{0) model is shown by the solid curve in 
Fig. 10. With the selection function given by the filter 
transmission curve, we found a best fitting value of ro = 
2.86 zb 0.33 h~^Mpc - shown by the dashed curve in Fig. 10. 
Given that these two estimates of ro are consistent with each 
other within the estimated uncertainties, it is clear that the 
selection function is not significant for the result. We note 
that the data-points appear to rise above the fitted model 
at small separations {0 < 0.15''), potentially indicating the 
beginning of the one halo term. Indeed the scale at which 
this turn-up in the correlation function is seen is compara¬ 
ble to that reported for z 3 LBGs by Hildebrandt et al. 
(2009). 

Such a low clustering measurement is consistent with 
previous measurements of the clustering of z ~ 3 

LAEs. Those shown in Fig. 10 have reported clustering 
lengths of ro = 2.51 q 7 Mpc (Gawiser et al. 2007) 
and ro = 1.70lo’46 Mpc (Ouchi et al. 2010), whilst 
Hayashino et al. (2004) did not present a measurement of 
the clustering length. As discussed, the Ouchi et al. (2010) 
sample is based on a comparable EW cut as our own sam¬ 
ple and so presents the most like-for-like comparison. Gom- 
pared to their "best" estimate of the clustering length, we 
measure a moderately higher ro, at the ^ 2 —3a level. Both 
measurements are subject to uncertainties in contamination 
and they report a "maximum" clustering length estimate 
of ro = 1 . 99 lo’ 55 - Similarly, our own measurement may 
be ’over-corrected’ for contamination based on the spectro¬ 
scopic non-detections, i.e. if the non-detections were in fact 
LAE too faint to be detected in the VIMOS observations. 
If this were the case, and we assume no contamination of 
our sample, then we would measure a clustering length of 
ro = 2.18 h~^ Mpc zb 0.26, somewhat more consistent with 
the previous measurement. 

The LAE clustering is low compared to both the normal 
star-forming and passive galaxy populations at z ~ 2 —3. In¬ 
deed, at z ~ 3, the LBG population of star-forming galaxies 
is found to have clustering lengths of ro ~ 4 Mpc (e.g. 
Adelberger et al. 2005a; Bielby et al. 2011, 2013). Given lit¬ 
erature results suggesting a strong link between galaxy stel¬ 
lar mass and median halo mass (e.g. Bielby et al. 2014; 
McGracken et al. 2015), such a low clustering length is 
strongly suggestive of LAEs being low-mass galaxies residing 
in low-mass dark matter halos. 


5.4 Dark Matter Halo Masses 

We now explicitly estimate median halo masses from our 
clustering results. Eollowing Bielby et al. (2014), we calcu¬ 
lated the mean masses of dark matter halos within the 
galaxy samples by using the clustering results. We used 
the formalism developed by Mo <V White (1996) which pro¬ 
vides a relationship between the halo-bias to the mean halo 
mass via Sheth et al. (2001). Based on the measured clus¬ 
tering length derived using the filter transmission profile, 
we find a clustering bias of 5 = 2.13 zbO.22 and halo mass 
Mdm = 
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Gawiser et al. (2007) reported a bias factor oi h — 
1.7 zb 0.4 and a median dark matter halo masses of Mdm = 
2 q 10 . 9 ±o .9 whilst Ouchi et al. (2010) reported a bias 

value oi b = 1.7zb0.8 and a halo mass Mdm ~ 6 . 7 lg^ 7 ° x 
10^° M©. We note that, although the cosmologies used in 
these papers and our own differ slightly, this only affects 
the halo mass values by a factor of ~ 10°'^. Ouchi et al. 
(2010) concluded that the average dark halo mass of LAEs is 
~ lO^i^i M© at z = 2 — 7 and our own measurement is more 
consistent with this assertion than their own results at z = 
3.1. As for the GALFORM model clustering result shown 
in Fig. 10, this corresponds to LAEs occupying dark matter 
halos with a median mass of Mdm ~ ]^q 11 - 3 ±o .3 Q^p_ 

proximately la higher than the measurement based on our 
observations. 

The median halo mass for LBGs is estimated to be 
~ 10^^ M©, about one order of magnitude larger than 
that of LAEs (Hamana et al. 2004; Ouchi et al. 2004b, 2005; 
Lee et al. 2006, 2009; McLure et al. 2010; Hildebrand! et al. 
2009). Indeed, from our own LBG measurements, we es¬ 
timate a median LBG halo mass at z ~ 3 of Mdm = 
10 ii. 6 ±o .2 ^ = 2.37 ± 0.21 (Bielby et al. 2013), 

larger than the LAE measurements at the ~ 2a level. 


5.5 LAE-LBG cross-correlation 

We now turn to the LAE-LBG cross-correlation, combin¬ 
ing the LAE candidates selected in this paper with the 
spectroscopically confirmed sample of LBGs presented by 
Bielby et al. (2013). Gombining the cross-correlation of the 
two samples with their respective auto-correlation samples, 
we may infer the relationship between the dark matter den¬ 
sity distributions that they each inhabit (see for example 
Ryan-Weber 2006; Tejos et al. 2014). 

From the Gauchy-Schwartz inequality (e.g. Gauchy 
1821; Mitrinovic et al. 2013), Adelberger et al. (2003) de¬ 
rived the following inequality for the auto- and cross¬ 
correlation functions of the two populations: 


^ab ^ ^aa^bb (9) 

where ^ab is the cross correlation and ^aa and ^bb are the 
respective auto-correlation functions. As discussed (in ref¬ 
erence to the galaxy-Hi cross-correlation) in Tejos et al. 
(2014), if the equality holds, the two samples are deemed 
to trace the same underlying mass distribution, and the rel¬ 
ative biases can be used to infer the ratio between the dark 
matter halo populations that the two populations occupy. 

We calculate the cross-correlation between the photo¬ 
metric LAE sample and a subset of the LBGs matching the 
criteria I^^lbg — 3.08| ^ 0.08 (giving 80 LBGs). Although 
this leaves a relatively small number of LBGs, it optimises 
the cross-correlation signal by excluding poorly correlated 
pairs at large separations along the line of sight. For the 
calculation, we use the Landy-Szalay estimator: 


w{6) 


{DaDbje)) - {DgRbiO)) - {RaPbiO)) + {RgRbiO)) 
(RaRbiO)) 


( 10 ) 

where DaDb is the LBG-LAE pair count, DaRb is the LBG- 
random LAE pair count, DaRb is the random LBG-LAE 
pair count and RaRb is the random LBG-random LAE pair 
count (all normalised by the number of randoms per galaxy). 
We use random catalogues with 20 x the number of randoms 
as galaxies in each case. 

The results are shown by the triangles in Fig. 11, 
alongside the LAE auto-correlation result (filled circles). 
We fit the cross-correlation result based on a single power- 



9 (arcmin) 

Figure 11. The LAE-LBG angular cross-correlation function cal¬ 
culated using the VLRS LAE and LBG samples (triangles). A fit 
to the data using a single power-law form (with a correlation 
length of ro = 3.29 ± 0.57) is shown by the solid curve. For com¬ 
parison, the LAE auto-correlation is reproduced (circles). 


law form in real-space, i.e. <f(r), transformed to w{0) based 
on the formalism outlined in Sec. 5.3 (but now incorpo¬ 
rating the two separate redshift distributions rather than 
just that of the LAE). The best-fitting result is shown 
by the solid curve in Fig. 11, corresponding to a correla¬ 
tion length of ro = 3.29 zb 0.57 /i“^Mpc (we do not fit 
for the slope, but take a fixed value of 7 = —1.8). This 
compares to a LAE auto-correlation clustering length of 
ro = 2.86 zb 0.33 /i“^Mpc and a LBG auto-correlation clus¬ 
tering length of ro = 3.83 zb 0.24 h~^Mpc (with a slope of 
7 = —1.6, Bielby et al. 2013). 

Evaluating Eq. 9 using the clustering amplitudes corre¬ 
sponding to r = 1 /i“^Mpc (i.e. r^'), we find ^ab/^aa^bb = 
1.28 zb 0.46, i.e. 1 . As discussed by Adelberger et al. 
(2005b); Tejos et al. (2014), this unity in the relation does 
indeed imply that the two photometric selections give sam¬ 
ples of galaxies that consistently trace the same underlying 
dark matter distribution and we can therefore use this result 
to provide a second constraint on the typical halo mass that 
our LAE sample occupies. Similar to the auto-correlation 
case, the correlation function can now be related to the un¬ 
derlying dark matter correlation function by: 


^LAE-LBG = 5lAe5lBG^DM (H) 

where ^lae-lbg is the LAE-LBG cross-correlation, 5lae is 
the LAE bias, 5lbg is the LBG bias, and ^dm is the un¬ 
derlying dark matter clustering as in the previous section. 
Based on this and our cross-clustering results, we calculate 
a further measure of the LAE bias of 5lae = 2.24zb 0.37. Us¬ 
ing this to estimate a halo mass as in the previous section, 
we find a value of Mdm = ]^q1i i±o. 4 xhe cross¬ 

correlation results corroborate the auto-correlation analysis 
with both dark matter halo estimates being consistent. 


6 DISCUSSION 

Our results, based on the largest z ^ 3 narrow band LAE 
sample yet surveyed and the first to combine multiple inde¬ 
pendent large fields, builds on the results presented by others 
(e.g. Gronwall et al. 2007; Gawiser et al. 2007; Ouchi et al. 
2008, 2010 ) to add to the overall picture of the nature of 
LAE galaxies at this epoch. 

At first sight, LAEs at z ~ 3 seem to be a high density 
population with a sky density of ~ 0.1 arcmin”^ in a redshift 
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range of width only Az ~ 0.05 to the various limits used in 
our survey. Taking the sky density of LBGs to be 1.8 
arcmin”^ in the 2.5 < z < 3.5 redshift range (Steidel et al. 
2003), this corresponds to a roughly similar sky density. On 
the other hand, comparing galaxy numbers at fixed contin¬ 
uum luminosity (e.g. Ouchi et al. 2008, Fig. 8 in this paper) 
we confirm that LAE are a relatively rare component of the 
galaxy population with abundances lOx lower than that 
of LBGs. 

Given their relative rareness, one might expect that 
LAEs would occupy rare high-mass dark matter halos based 
on a simple abundance matching approach. Our clustering 
results show a somewhat higher clustering signature than 
the previous best study of LAEs at this redshift (Ouchi et al. 
2010 ), however the clustering results still clearly show that 
the LAEs are relatively poorly clustered compared to av¬ 
erage star-forming galaxies at z 3 (e.g. Adelberger et al. 
2005a; Bielby et al. 2013; Durkalec et al. 2015). Indeed our 
measurement of LAE clustering implies a median halo mass 
of Mdm ~ 10^^ for which the number densities 

are predicted to be ~ 0.008 /i^Mpc“^ (integrating the mass 
function over the mass range Mpc) in a AGDM 

Universe (Jenkins et al. 2001; Murray et al. 2013). It there¬ 
fore follows (given our measured number density of LAEs 
of ~ 2 X 10“^ h^Mpc~^) that LAEs possess a relatively 
low occupation number (i.e. number of galaxies per halo) 
of << 1 and trace the dark matter distribution relatively 
sparsely (corroborating the results of Ouchi et al. 2010). It 
is common to link a galaxy sample’s properties as far as pos¬ 
sible to its host halo properties, but it is apparent that the 
galaxy properties that govern LAE selection (i.e. Lya equiv¬ 
alent width) are unlikely to be well correlated with the halo 
mass. Indeed, our results support the hypothesis that LAEs 
presence within dark matter halos is stochastic in nature 
(Nagamine et al. 2010; Ouchi et al. 2010). 

Our measured halo mass of Mdm = ]^q1i o±o-3 
is consistent with estimated LAE halo masses across a wide 
redshift range of 2 < z < 7 (i.e. the full range of redshifts 
probed by available LAE samples). Indeed, our measurement 
brings a tighter agreement with available measurements at 
other redshifts than the previous best measurement aX z ^ 3 
(i.e. Ouchi et al. 2010). As such, this enhances the assertion 
made by Ouchi et al. (2010) that LAEs exist as a stage in 
the evolution of at least some galaxies upon the host dark 
halo reaching a mass of Mdm ~ 10^^ 

The LAE sample as selected using the NB497 filter lies 
within the redshift range of LBGs selected based on their 
U — B {or U — G) colours. It is worth considering then how 
the galaxies identified by these two selections relate to each 
other. LBGs by definition are galaxies that possess a Lyman 
Break, however this is not intrinsic to the galaxies them¬ 
selves and so LAEs are also LBGs in the sense that they 
will have a Lyman break of some form, although it still may 
be too weak to satisfy the LBG cut (Gooke et al. 2014). Nev¬ 
ertheless, since the LAE number density at fixed continuum 
luminosity is ~ lOx lower than that of the LBGs, the LAE 
selection appears simply as an extreme subset of the LBG 
selected population. Indeed 40% of the LBG population has 
been shown to exhibit Lya in emission, with 25% of LBGs 
having equivalent widths of > 20 A (Shapley et al. 2003) 
and a tail in Lya equivalent width up to ~ 200 A. Glearly, 
at least some LAEs would fall into an LBG selected sam¬ 
ple and the numbers seem to approximately tally given that 
Gronwall et al. (2007) reported number densities of > 20 A 
equivalent width LAEs were approximately a third of the 
number densities of LBGs. Indeed, Eig. 8 of Shapley et al. 
(2003) shows that ~ 12% of their LBG sample have equiv¬ 
alent widths of > 60 A, which ties in very well with the 
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Figure 12. We show the clustering results as a function of limit¬ 
ing continuum magnitude for our LAE sample (large circles) and 
a range of 2 ; 3 LAE (circles) and LBG (triangles) samples. The 

dashed curve shows the predicted luminosity dependence based 
on a 1:1 relationship between luminosity and halo mass. 


number densities we find for our sample (as also noted by 
Ouchi et al. 2008 for their own 2 : ~ 3.1 EW > 60 A LAE 
sample). 

Our clustering analysis corroborates this picture of the 
LAE selected sample being an extreme subsample of the 
LBG population, whereby they show lower clustering due 
to our LAE sample probing a low-luminosity subset of the 
LBG population (at least compared to the mostly R < 2b 
LBG samples from which clustering measurements have 
been made). This is consistent with the UV luminosity de¬ 
pendence of the observed LBG clustering, whereby LBG 
samples selected with fainter UV luminosities are observed 
to have lower clustering signals (e.g. Giavalisco & Dickinson 
2001; Lee et al. 2006, 2009; Yoshida et al. 2008). In Eig. 12, 
we show the clustering results taken from a number of 
LBG (triangles) and LAE (circles) surveys including our 
own LBG clustering measurement and measurements of the 
clustering of our LAE sample when limited by a range 
of continuum magnitudes (filled black circles). The liter¬ 
ature z ~ 3 LAE results shown are from Gawiser et al. 
(2007); Ouchi et al. (2010), and the z ^ 3 LBG samples 
are from Giavalisco V Dickinson (2001); Adelberger et al. 
(2005a); Hildebrandt et al. (2009); Bielby et al. (2013). The 
dashed curve shows the predicted ro-magnitude dependence 
assuming a 1:1 ratio between UV continuum luminosity and 
halo mass, and using the Mo V White (1996) formalism to 
derive clustering lengths. 

We find a tentative continuum luminosity dependence 
in our LAE results, albeit only at the ^ la level. Taking our 
own LBG clustering results (which include the 5 fields used 
for the LAE study in this paper), our results are in broad 
agreement since our LBG limit is Rab ~ 24.8 and our lowest 
LAE continuum limit in Eig. 12 is Rab ~ 27 or a factor of 
8 in UV luminosity. Assuming that luminosity is approx¬ 
imately proportional to mass, given the above halo model 
this represents a factor of ~ 1.5 in bias. This compares to 
the observed factor of 1.35 ± 0.2 given the measured ro for 
LBG and LAE that we have measured above. The LAE clus¬ 
tering length of ro = 2.86 ± 0.33 /i“^Mpc is consistent with 
the decreasing trend of LBG clustering length measured by 
Hildebrandt et al. (2009). The VLRS LAE and LBG sam¬ 
ples show an Rum-'ro dependence consistent with this sim¬ 
ple model. At brighter magnitude cuts however, a steeper 
dependence is exhibited by the results of Adelberger et al. 
(2005a) and Hildebrandt et al. (2009). 

Larger LAE datasets are required to test if at the 
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Rab ~ 25 LBG limit, LAEs and LBGs have the same clus¬ 
tering amplitude, but at this point it seems that there is 
consistency with a model where LAEs can be treated as a 
subset of LBGs and where halo mass is proportional to UV 
luminosity in this combined population. The LBGs have a 
spatial number density that is consistent with their halo 
masses implied from clustering but the LAE number den¬ 
sity is too low, implying that they can only randomly sample 
their potential halo hosts. Erom the cross-correlation results 
and the equality of the ratio between the cross-correlation 
and the respective auto-correlations, it is also evident that 
although the fainter LAEs only random sample smaller mass 
halos, both the LBG and LAE populations trace the un¬ 
derlying dark matter distribution in a linear way but with 
significantly different bias. 


7 CONCLUSIONS 

In this paper, we have presented NB497 band photometric 
observations within 5 of our VLRS fields, taken using Subaru 
Suprime-Gam. Gombining these narrow band images with 
B, V and R deep imaging, we have made a photometric se¬ 
lection of z ~ 3.1 LAEs with an equivalent width limit of 

65 A. We have also made spectroscopic follow-up obser¬ 
vations of LAE candidates in 3 of the 5 fields, using VLT 
VIMOS. Two factors mean that this survey improves on pre¬ 
vious work on the study of LAEs at comparable epochs: 
firstly our choice of 5 independent survey fields minimises 
the potential effects of cosmic variance on our results com¬ 
pared to previous work; and secondly that our fields contain 
significant z ^ 3 LBG data has allowed us to perform a first 
cross-correlation analysis of the two populations. 

The major findings of our study are summarised as fol¬ 
lows: 

1. Our final LAE selection, based on the combination 
of B — ?nNB 497 with either V — mNB 497 or R — mNB 497 , gives 
a success rate of 78 ± 18% based on our spectroscopic ob¬ 
servations. The equivalent width limit of our final selection 
is ~ 65 A, whilst the flux limit is ~ 2 x 10“^^ erg cm^ s“^ 
(equivalent to L ^ 10^^ erg s“^). A catalogue of the photo¬ 
metric LAEs from the 5 VLRS fields presented here is avail¬ 
able at the Strasbourg Astronomical Data Genter (GDS, 
http://cds.u-strasbg.fr/ ). 

2. The spectroscopic observations produced 35 con¬ 
firmed LAEs at z ~ 3.1, with 23 of these falling within our 
final optimised selection criteria. We found no strong evi¬ 
dence for any contamination from low redshift On emitters 
(although < 5% of the detections remain ambiguous given 
the observational constraints). The only significant ‘contam¬ 
ination’ of the selection is from objects with no discernible 
emission in the observed spectra, pointing to either low Lya 
flux objects and/or interlopers (both likely entering the se¬ 
lection due to photometric uncertainties). 

3. We have calculated the NB497 number counts for 
our LAE sample (corrected for contamination) and find it 
to be consistent within the uncertainties with the previous 
observations of Gronwall et al. (2007, equivalent width limit 

20 A) and Ouchi et al. (2008, equivalent width limit 
65 A). Our selection is shown to probe the knee of the LAE 
Lya luminosity function based on an observed turnover in 
the narrow band number counts at ttinb ~ 24.5. 

4. We derived the R-band/rest-frame UV continuum lu¬ 
minosity functions of our sample of LAEs, again corrected 
for contamination. Our luminosity function lies at a level 

10 X lower than that of the z ^ 3 LBG luminosity func¬ 
tion, consistent with the LAE population being a small sub¬ 
set of the LBG population. At R < 26 (Muv ^ —19.6) our 


results are consistent with the fraction of strong emitters be¬ 
ing constant or marginally increasing towards fainter contin¬ 
uum luminosities, consistent with the results of Ouchi et al. 
(2008). Over a comparable magnitude range, we find our UV 
continuum luminosity function is significantly lower (by a 
factor of > 3) than that measured by Gronwall et al. (2007), 
however this appears consistent with their lower equivalent 
width threshold of ~ 20 A. 

5. We have measured the angular correlation function 

for our z ^ 3.1 LAE photometric sample in our 5 fields. Our 
results (corrected for contamination) are significantly lower 
than the observed clustering of LBGs at the same redshift. 
We measure a clustering length of ro = 2.86 ± 0.33 /i“^Mpc 
(assuming a slope of 7 = 1 . 8 ), which corresponds to a clus¬ 
tering bias of b = 2.13 ±0.22 and a median halo mass of 
2 q 1 i.o±o .3 This indicates that LAEs reside in low- 

mass dark matter halos, but given their number densities 
they have a low occupation number - i.e. only a small frac¬ 
tion of such halos actually host an LAE. 

6 . We measure the LBG-LAE cross-correlation function 
for the first time and find the results to be in agreement with 
the auto-correlation analysis, providing a useful consistency 
check on our results. Erom the cross-correlation analysis, we 
find a bias for the LAE population of 5 = 2.24 ± 0.37 and a 
typical dark matter halo mass of Mdm = ]^q1i i±o.4 

7. Thus our results for LAE number counts, clustering 
and cross-clustering are in agreement with a view that LAE, 
as selected here, are a low luminosity subset of LBGs that 
inhabit low mass dark matter haloes with a small occupancy. 
Gomparison to current semi-analytical models show good 
agreement in the narrow band counts, but poorer agreement 
in the LAE clustering properties. It will be interesting to 
see how future changes to semi-analytical models can better 
match the observational results. 
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APPENDIX A: VLT VIMOS ID AND 2D 
SPECTRA OF CONFIRMED LAES 

The figures show the observed-frame spectra taken using 
VLT VIMOS with the HR_Blue grism as part of our spec¬ 
troscopic follow-up observations. All the spectra shown have 
been identified as LAEs with significant emission within the 
NB497 filter passband. 
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Figure Al. VIMOS HR_red spectra (black solid curve) for confirmed LAEs in the QSO B2359+0653 field. The grey filled region shows 
the transmission profile of the NB497 filter. In addition, the insets show the corresponding 2D spectra showing only the wavelength 
region covered by the filter. The dashed lines show the noise estimated from the standard deviation of the background signal across the 
slit. 
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Figure A2. As for Fig. Al, but for QSO J0124+0044 field. 
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Figure A3. As for Fig. Al, but for LBQS Q0301—0035 field. 
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